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ABSTRACT 

> 

00 \ We investigate techniques that can be used to determine ages of starburst regions 

^ ■ containing populations beyond their early nebular phase. In particular, we study the 

^ '. strength of the Call triplet (AA 8498, 8542, 8662 A) and the CO index (2.31-2.40 //m 

O ! band) using synthetic models as the starburst evolves. For an instantaneous burst of 

star formation both of these absorption features remain strongest between 7-14 Myr 
corresponding to the red supergiant population. The detailed evolutionary behavior 
p • of the starburst is strongly metallicity dependent. Low metallicity starburst models 

successfully reproduce the distribution of equivalent widths of Ca ll triplet with age in 



o 

H \ Large Magellanic Cloud clusters. The clusters in the red supergiant phase strongly favor 

the stellar evolutionary models incorporating mass-loss rates higher than the standard 
values. We suggest usage of diagrams involving Ca ll triplet equivalent width, CO index 
and nebular recombination lines to infer the history as well as age of starburst regions. 
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1. Introduction 

Detection of Call triplet (CaT) lines in a cir- 
cumnTiclcar Hll region of the starburst galaxy 
NGC3310 by Tcrlcvich ct al. (1990) offered a 
new tool to study aging starburst systems. The 
CaT lines originate in the atmospheres of cool 
stars with the absorption being the strongest in 
red supergiants (RSGs). It is known from the 
stellar evolutionary models at solar metallicity 
(Z0=O.O2) that O stars in the mass range 15- 
25M0 evolve into RSGs in 7-14 Myr, and hence 
the detection of the CaT lines in starburst sys- 
tems is often used to infer the presence of a stellar 
population of around 10 Myr age. However, the 
CaT line equivalent widths of cool giants are non- 
negligible and hence the detection of the triplet 
lines is not an unambiguous signature of the pres- 
ence of RSGs. 

Bica, Alloin &: Santos (1990) provided a so- 
lution to this problem by obtaining the CaT 
equivalent widths of blue clusters with well de- 
termined turn-off ages in the Large Magellanic 
Cloud (LMC). They found that only clusters 
with an age '^10 Myr have CaT equivalent widths 
greater than 6A; the older clusters having values 
around 4A. However, their observed relation has 
only limited use in estimating the age of starburst 
regions due to the fact that the evolutionary his- 
tory of massive stars, and hence the presence of 
RSGs, is a strong function of metallicity (Maeder 
1991), and that starbursts are found in a wide 
range of metallic environments. Thus a quan- 
titative interpretation of the observed equivalent 
widths in starbursts requires studies of the evolu- 
tionary behavior of starbursts at different metal- 
licities. We carry out such a study here with 
special emphasis on reproducing the observed re- 
lations of Bica et al. (1990). We also study the 
evolutionary behavior of the CO index, which is 
also sensitive to the population of cool stars, to 
investigate its use as an age indicator in conjunc- 
tion with the CaT equivalent widths. As this is 
one of the first attempts in synthesizing these 



spectral features (see Leitherer et al. 1996 for 
a list of on-going attempts) we mainly focus on 
understanding the basic evolutionary properties 
and its metallicity dependence. The synthesis 
technique and the stellar data used are described 
in Sec. 2. Evolutionary results are presented in 
Sec. 3 and are discussed from the context of ob- 
servational data of Bica et al. (1990) for LMC 
clusters in Sec. 4. The results are discussed from 
the point of view of age estimations of starburst 
regions in Sec. 5. 

2. Synthesis model and Stellar Data 

We use the conventional evolutionary popu- 
lation synthesis technique to compute the star- 
burst related quantities. The synthesis code we 
use in the computation is described in detail in 
Mayya (1995) and the results are compared with 
other existing codes in Chariot (1996). The ma- 
jor change in the present code is the implementa- 
tion of the isochrone synthesis technique (Chariot 
&; Bruzual 1991) which has the effect of reduc- 
ing the noise on the computed quantities. The 
equivalent width of the CaT absorption lines is 
computed for an instantaneous burst (IB) of star 
formation using the formula, 

I Lca,T{m;t; Z)!^,'" dm 

EWcaT(imf,i,Z) = 

/ Lsmoim;t; Z)m-°' dm 

mi 

where LcaX is the strength of the CaT absorp- 
tion feature and Lggoo is the underlying contin- 
uum luminosity per unit A at the rest wavelength 
of 86OOA. The initial mass function (IMF) slope 
a = 2.35 (Salpeter 1955) is chosen with mass 
limits mi = 1 and m„ = 100 in solar units. 
The CO index is computed by a similar formula 
with Leo replacing LcaT and the corresponding 
continuum wavelength being 2.36^m. It is ex- 
pressed in magnitude units following the conven- 
tional definition. Stellar evolutionary models are 
used to obtain the spectral type and luminosity 
class for a star of mass m at evolutionary stage t 
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and metallicity Z . Theoretical model spectra are 
then assigned to each of these spectral types to 
enable us to compute the continuum luminosities 
at 86OOA and 2.36//m. The spectra from Kurucz 
(1992) are used for effective temperatures of stars 
above 3500 K. At lower temperatures we use the 
library of cool star spectra compiled recently by 
Lejeune, Cuisinier h Buser (1996). 

We use stellar evolutionary models from Geneva 
for starbiirst evolutions at mctallicitics Z = 
0.001, 0.008, 0.02 and 0.04. These metanici- 
ties encompass a wide range of starburst systems 
from the metal poor blue compact galaxies to 
the metal rich nuclear starbursts. Geneva mod- 
els referred as standard (Schaller et al. 1992; 
Schaerer et al. 1993a,b) assume moderate over- 
shooting of the convective core, and use empiri- 
cal relations for expressing mass loss rate (M) for 
different masses at different metallicities. There 
are known to be difficulties in explaining all the 
observables with one set of assumptions (Maeder 
1994; Langer & Maeder 1995). Since our interest 
here is to compute quantities related to RSGs, 
we looked for those tracks which correctly repro- 
duce the observed ratios of blue to red super- 
giants (B/R). Meynet (1992) has shown that the 
standard tracks reproduce the observed B/R at 
solar metallicity, while at lower metallicity, mod- 
els with mass loss rates enhanced by a factor 2.2 
reproduce the observed B/R better. Hence we 
used standard tracks at Z >Zq and Meynet et 
al. (1994) models with enhanced mass loss rates 
(by a factor of two over standard M) at lower 
metallicities. 

2.1. Stellar data for Can triplet and CO 
index 

The most comprehensive investigation on the 
dependence of the equivalent width of the CaT 
with spectral type, luminosity class and metal- 
licity of stars was carried out by Diaz, Terlevich 
& Terlevich (1989) (DTT hence forth). They 
defined the CaT equivalent width, EW(CaT), 
as the sum of the equivalent widths of two of 



the brightest lines (AA 8542, 8662 A) and pro- 
vided an expression relating EW(CaT) with log 5 
and log(Z/Z0). For metal rich stars {Z >Zq) 
the EW(CaT) increases with decreasing surface 
gravity and for supergiants, it increases with in- 
creasing metallicity. The dependence on effective 
temperature was found to be negligible for their 
sample of cool stars. Independent observational 
studies by Zhou (1991) and Mallik (1994) sup- 
port the above findings. Theoretical computa- 
tions of Jorgensen, Carlsson & Johnson (1992) 
(JCJ henceforth) successfully reproduced all the 
basic trends established by DTT. However they 
found the theoretically based relations involving 
higher order terms in logg are better fits to the 
observed data than the linear relation of DTT, 
which underestimates the EW(CaT) by ~ 2-5A 
in metal rich red supergiants. Hence we use JCJ 
relation for Z >Zq and DTT relation at lower 
metallicities. 

Data on the CO band is not as complete as 
that for the CaT. The best published data-set so 
far is that of Kleinmann & Hall (1986) for a sam- 
ple of 26 stars covering a wide range of spectral 
type and luminosity classes, but at solar metallic- 
ity only. Doyon, Joseph & Wright (1994) defined 
the CO index based on the spectra in the above 
sample and obtained relations for the CO index 
as a function of effective temperature for dwarfs, 
giants and supergiants separately. These make 
up the basis for our computations. The above 
defined index (COgp) is related to the photomet- 
rically measured index (COph; Frogel et al. 1978) 
by COsp =1.46 COph - 0.02. 

3. Synthetic results for Can triplet and 
CO index 

The evolutionary behavior of the EW(CaT) at 
the four metallicities is shown in Fig. la. The fol- 
lowing important stages can be recognized during 
the evolution of solar metallicity starbursts. 
Stage 1. EW(CaT) rises to detectable levels 
in starbursts, 5 Myr after the onset of the IB, 
reaching the highest values of ~llA between 7- 
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14 Myr. 

Stage 2. The EW abruptly drops to around 5 A 
at ~15 Myr, remaining at these low values up to 
around 20 Myr. 

Stage 3. The EW rises again above 6.5A at 
25 Myr, showing a secondary peak at 60 Myr. 
Stage 4. The EW attains an asymptotic value of 
around SA after 100 Myr. 



IB model with IMF: Mi=l. M„=100, a=2.35 
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Fig. 1. — Evolution of (a) the equivalent width of 
Call triplet and (b) the CO index for an instan- 
taneous burst of star formation for the assumed 
Salpeter initial mass function, at the four chosen 
metallicities. Supergiant phases of massive stars 
at 2Zq (top) and Zq metallicities are shown by 
the parallel horizontal bars, with the numbers in 
the middle denoting the stellar masses in Mq. 

The evolutionary behavior at Z =Zq/2.5 is 
qualitatively similar to Zq models. However 
the peak values, as well as the contrast between 
Stages 1,3 and 4 above, are lower. The trend of 
decreasing contrast continues at lower metallicity 
resulting in steady values of around 2A through- 



out the Stages 3 and 4 for Zq/20 model. The 
most significant feature during the evolution of 
low metallicity starbursts is the absence of the 
primary peak at 10 Myr (Stage 1). At intermedi- 
ate metallicities (Zq/5; not plotted), the 10-Myr 
peak is present (2A) but at values lower than at 
Stage 3 (5A). At higher than solar metallicities. 
Stage 2 is missing resulting in a smooth decrease 
of the EW from Stage 1 to the end of Stage 3. 

The Stages 1-4 as well as their dependence 
on metallicity can be understood by identifying 
the individual phases, which contribute signifi- 
cantly to the EW(CaT), during stellar evolution. 
The cool giant and supergiant phases of different 
stellar masses are depicted in Fig. 1 by horizon- 
tal bars at the upper part of the diagram. The 
primary peak at ~ 10 Myr at all the metallicities 
is directly related to the cool RSG population of 
progenitor masses between 15-25Mq. The su- 
pergiants at the lowest metallicity we have con- 
sidered, spend most of their post-main sequence 
evolution as blue supergiants, and hence fail to 
show the peak at 10 Myr. 

Note that the horizontal bars depicting the 
cool giant and supergiant phases have a discon- 
tinuity for stars of mass < I2M0 at solar metal- 
licity. This break is due to the qualitatively 
different post-main sequence evolution of 12Mq 
and lower mass stars for Z <Zq. These stars 
go through the well-established "blue loop" in 
the Hertzsprung-Russel diagram during the core 
helium burning phase, returning to the redder 
phase only at the end of core helium burning. 
Thus, there is a sudden reduction in the num- 
ber of RSGs in the starburst, which explains the 
Stage 2. "Blue loops" during core helium burn- 
ing is a common property of stellar evolutionary 
models for Z <Zq, and hence Zq/2.5 model be- 
haves qualitatively similar to the Zq model. On 
the other hand, "blue loops" are absent for higher 
metallicity models and hence there is a smooth 
transition from the RSG-dominated phase to the 
red giant and asymptotic giant branch domi- 
nated phases at Z = 2Zq. 
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Successively lower mass stars spend a lesser 
fraction of their post-main sequence lifetime in 
"blue loops" , resulting in the accumulation of red 
giant and asymptotic giant branch stars, leading 
to the secondary peak at 60 Myr for models. 
The adopted scaling between mass-loss rate and 
metallicity and also the corresponding phases for 
a star of given mass being warmer at lower metal- 
licities result in the net decrease of the number 
of stars contributing to CaT absorption at lower 
metallicities. This is the reason for lower peak 
values as well as lesser contrast between differ- 
ent Stages at lower metallicities. After 100 Myr, 
the EW(CaT) is controlled by the cool giant and 
dwarf stars of mass < 5Mq. 

Evolution of the COgp index is shown in Fig. lb. 
Note that the metallicity dependence in this plot 
arises only due to the dependence of stellar evo- 
lution on metallicity; there is no data-set on its 
dependence on line formation. The evolution- 
ary behavior, as well as its dependence on metal- 
licity, can be understood by the same physical 
phenomenon as that for the EW(CaT). However 
the following quantitative differences in the four 
Stages are noteworthy. 

Stage 1. The RSG peak is narrower for Z <7iq 
(8-12 Myr). 

Stage 2. The CO index rises faster after the drop 
following the "blue loop" phase. 
Stage 3. The secondary maximum occurs at 20- 
30 Myr, decreasing smoothly then onwards up to 
around 100 Myr. 

Stage 4. The CO index do not reach asymptotic 
values beyond 100 Myr. 

These differences arise firstly, due to a stronger 
dependence of the CO index on the effective tem- 
perature and secondly, due to the same stars con- 
tributing to the continuum as well as the CO 
absorption strength. The CO index of the star- 
burst is essentially that of the most luminous cool 
star, as these cool stars contribute both to the 
line depth as well as the continuum. In contrast 
the continuum at CaT has a significant contribu- 
tion from hotter stars and hence the EW(CaT) 



evolves differently than that of the CO index. 

4. Comparison v^'itli LMC clusters 

The blue clusters of the LMC provide a unique 
opportunity of checking the evolutionary results 
discussed above, before these results can be used 
to derive ages of starburst systems. The spe- 
cial advantage of these clusters is the availability 
of well-determined turn-off ages. As discussed 
earlier, Bica et al. (1990) had obtained the 
EW(CaT) for each cluster, treating each clus- 
ter as a single entity. Bica et al. also com- 
piled the CO indices and broad band colors for 
each of these clusters, which are ideal for com- 
parison with our results. We superimpose our 
model results on these observed data in Fig. 2. 
The strength of the fainter A8498A line has been 
subtracted from the Bica et al. value to covert 
the observed values to the more commonly used 
definition of DTT. Photometrically derived CO 
indices are converted into COgp using the rela- 
tion given in sec. 2.1. The available J — K and 
B — V colors are also plotted in the last two pan- 
els. The solid and dashed lines correspond to 
models using enhanced and standard mass loss 
rates respectively, both at Z0/2.5 metallicity. 

The ages of the clusters with the highest 
EW(CaT), CO index and J — K color agree very 
well with the RSG dominated bump (Stage 1) 
in our model. The absolute values of the RSG 
peak with stellar evolution incorporating stan- 
dard mass loss rate is too low to account for the 
observations, confirming the under-estimation of 
RSGs in these models. Usage of enhanced mass 
loss rate models correctly reproduces the EW(CaT) 
and J — K color. The observed CO index is 
still higher by around 0.08 mag, which might 
be caused by errors introduced in converting 
the photometrically observed indices to spectro- 
scopic indices. The larger scatter in this plot 
compared to the other plots also suggests larger 
errors in the observed CO indices. Thus we con- 
clude that the stellar evolution with enhanced 
mass-loss rates, and not the standard, is a better 
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ZO/2.5 IB model with IMF: M,=l, M„=100, a=2.35 
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Fig. 2. — Comparison of the computed quanti- 
ties at Z =Zq/2.5 with the observational data of 
Bica et al. (1990) for the LMC clusters. Stel- 
lar evolutionary models incorporating twice the 
standard mass-loss rate (solid line) reproduce the 
observed characteristics better. 

representation of the observed values in LMC. 

Older clusters have EW(CaT), CO index and 
J — K as predicted for Stage 3 in our model. 
However only the lower boundary of the observed 
B — V colors agree with the model predictions for 
Stage 3 with the majority of the clusters having 
colors around 0.1 magnitude redder. A similar 
trend is also seen in U — B colors (not plotted), 
which is most likely due to the internal reddening 
in the clusters. It is interesting to note that the 
observed B — V (and also U — B) colors during 
the RSG phase are bluer than the predictions of 
the enhanced mass loss rate models. 

5. Discussions and Outlook 

Having studied the evolutionary behavior of 
the CaT equivalent width and the CO index, we 
now discuss their usage in estimation of ages of 
starburst regions from observed data. Diagrams 
involving dimensionless quantities play an impor- 
tant role in this direction. Most of such attempts 
hitherto arc based on the presence of massive O 
stars or Wolf-Rayet stars and hence are applica- 



ble only as long as the starburst nebula is trace- 
able. The end of nebular phase coincides with 
the beginning of the RSG phase for an IB, as the 
ionization is mostly controlled by stars having 
masses higher than the progenitors of the RSGs 
(~ 3OM0). This is illustrated in the top panel 
of Fig. 3. It can be seen that by the time the 
EW(CaT) attains the peak value the Ha emis- 
sion equivalent width has dropped below the de- 
tection limit. We investigated various diagrams 
involving EW(CaT), CO index and colors for the 
four metallicities, and found the EW(CaT)-CO 
index plane to be the best suited for the purpose 
of age determination. The loci of IB and contin- 
uous star formation (CSF) models aX Z = TLq 
and Zq are shown in Fig. 3. Selected ages in Myr 
are indicated along the locus in this plot. 

In the EW(CaT) vs EW(Ha) plot, the IB 
and CSF models follow quite different tracks 
and hence can be most effectively used to in- 
fer the history of star formation; e.g. regions 
with detectable EW(CaT) and emission EW(Ha;) 
greater than 300A can be described only by CSF. 
For CSF models asymptotic values are reached 
in ~15 Myr. Once the history of star forma- 
tion is known, the EW(CaT) vs CO index plot 
can be used to derive accurate ages of starburst 
regions. The ambiguity in the estimated age in- 
creases with decreasing metallicity for IB models, 
with the result that similar plots for Z <'Lq have 
limited use. However our models give a definite 
upper limit for low metallicity systems. The blue 
compact galaxies with Z <Zq/^ are not expected 
to have EW(CaT) and CO index values above 5A 
and 0.1 dex respectively. These values are insen- 
sitive to mass limits as long as > 3OM0 and 
mi < 3M0, and are only weakly sensitive to the 
IMF slope. 

Before we end, we point out the limitations 
of the present work and discuss where further 
improvements are necessary. Our main purpose 
in this work was to investigate the important 
stages during evolution of EW(CaT) and CO in- 
dex in starburst regions and their metallicity de- 
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IMF: Mi=l, M„=100, a=2.35 
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Fig. 3. — Dimensionless diagrams involving ob- 
servational quantities for Z =Zq and Z =2Zq 
models with instantaneous burst and continuous 

star formation (CSF) at a constant rate. The 
numbers next to tick marks are starburst ages 
in Myr. For the CSF model values stabilize at 
~15 Myr. 

pendence, and accordingly we kept the starburst 
model at the simplest level. Observed values in 
nuclear starbursts may have some contribution 
from the bulge light (see e.g. Garcia- Vargas et 
al. 1997), and hence usage of Fig. 3 for nuclear 
starbursts may require a separate treatment of 
the bulge. However Fig. 3 is directly applicable 
for high metallicity giant Hll regions in disks of 
galaxies. Our next step is to study the dilution 
of the computed quantities function of frac- 
tional contribution of the bulge light, which will 
enable direct comparison with observed data of 
starburst nuclei. 
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